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Abstract 

We provide a systematic and thorough exploration of lepton flavor models in which 
the solar mixing angle is related to the golden ratio. For scenarios in which the solar 
mixing angle is given by the inverse cotangent of the golden ratio, we demonstrate that 
A§ is the smallest non-Abelian finite group that contains all of the symmetries necessary 
to enforce this specific lepton mixing pattern. Within this context, we propose two 
lepton flavor models that yield this mixing pattern through the breaking of A§ at 
leading order to the Klein four subgroup in the neutrino sector. Both models have 
triplet embeddings of the lepton doublets as well as the charged lepton singlets. In 
the charged lepton sector, the residual symmetry is Z§ in the first model, while in 
the second model, A§ is broken completely at leading order. For the second model, 
the reactor mixing angle vanishes at leading order and is of the order of the square of 
the Cabibbo angle at next-to-leading order, which is allowed by the global analysis of 
current lepton data. 



1 Introduction 



Neutrino detection experiments have now provided definitive evidence that the solar and 
atmospheric neutrino anomalies can be explained by neutrino oscillations. Global fits of the 
data for the mass-squared differences and the mixing angles of the Maki-Nakag aw a- S akat a- 
Pontecorvo (MNSP) matrix p]j2] assuming three families of neutrinos have been provided by 
several groups [3H7] ; we show the results of the latest two independent global fits in Table [TJ 
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Table 1: Summary of neutrino oscillation parameters from two global data fits [2l[5|. The upper (lower) row 
corresponds to normal (inverted) neutrino mass hierarchy in [5] , and the results from the new reactor fluxes 
are shown in parentheses for the fit of [S]. 



The neutrino oscillation data have revealed great differences between the quark and 
lepton mixing angle patterns. The atmospheric angle 023 is compatible with maximal mixing, 
but the error bars admit relatively large deviations. The solar angle 0i2 is more precisely 
measured and is known to be large but not maximal at the 5a level. The reactor angle 0i3 
has recently been the source of great excitement since T2K [8] and MINOS [9] each reported 
evidence for a relatively large 0i3 at the level of 2.5a and 1.7a, respectively. The first result 
from Double Chooz experiment also implies a non-zero 0i3 [TU] . The global fits quoted above 
that include these new results indicate 0i3 > at about the 3a level, but the magnitude of 
013 still suffers from large uncertainties. Future measurements of 0i3 at Double Chooz [11] 
and Daya Bay [12] will provide important clues in the search for a greater understanding of 
the fermion mass puzzle of the Standard Model (SM). 

In contrast to the small quark mixing angles, the large lepton mixing angles suggest 
a model-building paradigm based on discrete non-Abelian family symmetries (for reviews, 
see [13]). Lepton models have been constructed based on the tetrahedral group T (which is 
isomorphic to A4) [T3], the binary tetrahedral group T' [15J, A(3n 2 ) and A(6n 2 ) [16J, the 
semidirect product of Z 3 and Z 7 [XT] . Z 2 x Z 2 [H], S4 [19], S3 [20], the semidirect product 
of At and S3 [21], PSL{2, 7) [22], the quaternionic symmetries [23], the dihedral symmetries 
D n [2UI25], ^13 [2S], and the icosahedral group / (which is isomorphic to A 5 ) J2ZH2H] (an 
extension to I' for the quark sector can be found in [30]). Many scenarios lead to the 
well-known Harrison-Perkins-Scott [HPS] "tri-bimaximal" mixing pattern [31], in which at 
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leading order #23 is maximal, #13 is zero, and the solar angle is given by cot 0i2 = V2. We 
note that the tri-bimaximal matrix appeared earlier in the context of quark masses [32] . 
However, there are also alternative mixing patterns that are also in agreement with current 
data, but differ according to their leading order predictions for the solar mixing angle #12. 
Recently, a number of intriguing scenarios have also been proposed in which # 13 ^ at 
leading order [33H39] . in response to the recent possible hints of nonvanishing #13. 

In this paper, we focus on scenarios in which at leading order #23 is maximal, #13 = 0, 
and #12 is governed by the golden ratio, (j) g — (1 + v5)/2. In this context, there are two ideas 
that have been proposed. The first, which we call "GR1" mixing, is the hypothesis that 

COt 6l2 = (j)g, (1) 

which was first proposed in jlQ] and later explored in detail within a Z 2 x Z 2 flavor symmetry 
in [TS]. In this case, sin 12 = y/h/{h<t> g ) ~ 0.276, which is very close to the 3a lower boundjj 
As first suggested in [TS] , the icosahedral symmetry group I ~ A 5 is a natural candidate as 
a flavor symmetry that can yield this mixing pattern, since the golden ratio characterizes 
several geometrical properties of the icosahedron. As a result, A§ models based on this idea 
have been proposed in [27H29] . The second, which we call "GR2" mixing, is the relation first 
suggested by jH], in which 

cos£ 12 = ^. (2) 

Hence, in this case sin 2 # 12 = v / 5/(40 9 ) ~ 0.345, which is close to the 2a upper limit. A 
model realization of this mixing pattern based on dihedral symmetry groups was given in [25] . 

The purpose of our work is to provide a systematic exploration of the golden ratio mixing 
hypotheses. For both of the GR1 and GR2 mixing patterns, we examine the question of what 
is the minimal group that naturally contains all of the symmetries of the lepton matrices 
within each scenario. In the case of GR1 mixing, we find after a systematic analysis that the 
minimal group is indeed A 5 . In the case of GR2 mixing, there is no discrete group (up to a 
very large order) that contains the complete set of symmetries of the associated lepton mass 
matrices. Focusing our attention on the case of GR1 mixing, we then provide a comprehensive 
exploration of lepton flavor model building based on A§, with the key input that the lepton 
doublets and charged lepton singlets of the SM are embedded in triplet representations of 
A 5 as opposed to singlet embeddings. Within the framework of supersymmetric models, we 
obtain two models of this type that differ in terms of their respective residual symmetries in 
the charged lepton sector. In one model, A 5 is broken to Z 5 in the charged lepton sector and 
to Klein four in the neutrino sector, resulting naturally in GR1 mixing. However, fine-tuning 
is needed to obtain the charged lepton mass hierarchy. The second model improves on this 
scenario by breaking this residual symmetry. 

This paper is organized as follows. In Section 2, we present our investigation of the 
minimal family symmetry group required to produce the GR mixing patterns. As this 
procedure yields the minimal group of A 5 for GR1 mixing and no obvious solution for GR2 
mixing, in Section 3 we present basic group theoretic aspects of A 5 that are useful for flavor 
model-building, including the equivalency classes and the explicit subgroups. In Section 4, 

1 sin 2 0i2 ~ 0.276 is slightly above the 3cr lower bound of the first global fit, while it is around the la 
lower limit of the second fit. 
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we turn to A§ flavor model-building in the lepton sector. We summarize and conclude in 
Section 5. In the Appendix, further practical details of the group theory of A 5 are given, 
including the explicit representation matrices and the Clebsch-Gordan coefficients. 



2 Minimal Flavor Symmetry for Golden Ratio Mixing 

In this section, we address the question of finding the minimal discrete non-Abelian 
flavor symmetry group that can naturally realize the golden ratio mixing hypotheses. In this 
context, the minimal group is defined as the smallest finite group that contains as elements all 
of the symmetries of the neutral and charged lepton mass matrices. Given these symmetries, 
we then solve the inverse of the usual issue of finding group representations: given specific 
representation matrices, it is necessary determine the group such that these matrices form 
a faithful representation of the group (i.e., one in which the kernel of the representation is 
just the trivial subgroup consisting of the group's identity element). 

The first step is to determine the symmetries of the mass matrices. We start with the 
charged lepton sector. We choose to work in a basis which fhg = m £ mi is diagonal, where 
mn is the charged lepton mass matrix. Hence, it is straightforward to see that the symmetry 
transformation matrix Gg, which is determined by the condition G\rhgGi = frig, is a diagonal 
and non-degenerate 3x3 unitary matrix (i.e., a diagonal and non-degenerate phase matrix). 
This ensures that rhg remains diagonal since it has non-degenerate eigenvalues. 

Turning to the neutrino sector, with this choice of basis the neutrino mass matrix can be 
reconstructed in terms of its complex mass eigenvalues mi >2i3 as follows: 

m„ = U*dia,g(m 1 , m 2 , m 3 )W . (3) 

Here U is the MNSP matrix (in a specific phase convention) with 6 12 unspecified, # 23 = vr/4, 
and #13 = 0: 

/ C12 S12 \ 

U = - 8i2 /V2 c 12 /V2 1/V2 , (4) 
V -S12/V2 c 12 /V2 -1/V2 J 
in which s 12 = sin #12 and c 12 = cos#i 2 H The neutrino mass matrix then takes the form 



m. 



/ m r 2 _l_ m c 2 m 1 -m 2 m 1 ~m 2 \ 

I mic 12 + m 2 s 12 Si 2 Ci 2 s i2 Ci 2 ^ 

- zza ^f M Si 2 Ci2 \m\s\ 2 + \m 2 c\ 2 + |m 3 \m x s\ 2 + \m 2 c\ 2 - |m 3 

\ ~ m 7T lsi2Cl2 \ m A + \ m 2c\ 2 - |m 3 \mxs\ 2 + \m 2 c\ 2 + |m 3 j 



(5) 



We turn now to the GRl mixing pattern, in which cot 9\ 2 = <p g (a similar analysis can 
be found in [2H]). In this case, the neutrino mass matrix is given by 

m GR1 

b 9 $9 

(6) 





2 We absorb the Majorana phases in the mass eigenvalues mi, rather than in the MNSP matrix. 
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The 2 — 3 symmetry of Eq. fl6]) is due to the maximal atmospheric angle and the vanishing 
reactor angle. In addition, the following relation is satisfied: 



GR^ 



111 



V2\ 



m 



GRU 



112 



K nf \ u ) 2 2 + (m u ) 2 3- 



(7) 



To obtain the symmetry properties of m„ , one seeks the unitary transformations Gi that 
satisfy the condition Gjm^ R1 Gi = m^ R1 . After a straightforward computation, we find (in 
addition to the identity transformation) three unitary transformation matrices, G\ 23H 




Go 




V2 
-lid) 



s/2 



-1/ 



which satisfy the relations 

Cr? = 



GiGj 



GjGi = Gk with i 7^ j ^ k, 



(8) 



(9) 



Therefore, the symmetry group of the neutrino mass matrix m^ R1 is the Klein four group 
K 4 ~ Z 2 x Z 2 (the result that the symmetry group is Z 2 x Z 2 was also found in [29]). 

For the GR2 mixing pattern with cos #12 = (/> g /2, an equivalent reconstruction of the 
neutrino mass matrix results in 



m 



GR2 




3 - 
3 - 



m 2 1 

This matrix, which is also clearly 2 — 3 symmetric, satisfies the following sum rule: 





(0 







)+?! 





1 






^ 


-1 





(10) 



(m^ 2 )n + 2 




OR2>, 



112 



r _GR2N 1 / GR2' 
J 22 + K , 



23- 



Following the same procedure as above for determining the symmetry transformations G\ 




m 



GR2 



we find the following set of three unitary matrices: 



^-6 

2 <rg 



2 



VE 1 

2 </>„ 



%/5 1 
2 </>, 



/ 



f0 9 



2 ^3 
y 7 ! 1 



VE 1 

2 <A 



(12) 



3 We have chosen a convention for the overall sign of each Gi such that its determinant is equal to unity. 
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The G\ satisfy the relations given in Eq. (jUJ), so the symmetry of the neutrino mass matrix 
m^ R2 is also the Klein four group. Note also that G' 3 = G 3 ; this is as expected since it is 
determined solely by the requirements that #13 = and 6*23 = 7r/4. 

Based on these results, we apply the following procedure to find the minimal discrete 
family symmetry group for these mixing patterns. Denoting this underlying family symmetry 
group at high energies as Q, the above considerations tell us that Ge and Gi j2j3 {G[ 2 3 ) should 
be elements of this group for GR1 (GR2) mixing in order that they are preserved in the theory 
at lower scales at leading order (LO). For the case of finite groups, this implies that there 
must be some integers n and rrii such that 

G n t = {GiG t ) mi = 1, n > 3, (13) 

where the restriction on n results from the requirement that Gt is nondegenerate, and an 
identical relation holds for the primed version G' 123 . The group generated in this way for 
each Gi must be a subgroup of the full symmetry group Q. Hence, the procedure is first to 
determine possible solutions for Gi that result in the relation of Eq. (fl3|) with integer values 
for n and m 8 given Gi or G[, and then to examine their possible embeddings in discrete 
groups. Working in the diagonal basis for the charged leptons leads to no loss of generality 
because the results obtained in different bases are related via similarity transformations. 
Furthermore, the 2 — 3 symmetry of the neutrino mass matrices indicates that identical 
results are obtained for choices of Ge with the (22) and (33) elements of Gt interchanged. 

Let us consider the case of GR1 mixing as the first example. An explicit computation 
yields the result that for n = 3 (G 3 = 1) or n = 4 (G 4 = 1), the resulting family symmetry Q 
must be an infinite group, since there is no integer m for which (GiGe) m = 1 or (G 2 Ge) m = 1. 
However, n = 5 (Gf = 1) yields ten nontrivial solutions, which are tabulated below: 

• G t = (l,p,p 4 ): (G X G,) 3 = 1, (G 2 G e ) 5 = 1, (G 3 G e ) 2 = 1. 

• G £ = (l,p 2 ,p 3 ): (G^ef = 1, {G 2 G,f = 1, (G 3 G e ) 2 = 1. 

• G e = (p,l,p 2 ): (G 1 G i ) 15 = 1, (G 2 Ge) 5 = 1, (G3GV) 10 = 1. 

• G t = (p 4 , l,p 3 ): (GiG e ) 15 = 1, (G 2 G e f = 1, (G3GV) 10 = 1. 

• Gt = (p 2 ,P,p 3 ): (G1GV) 15 = 1, (G 2 G e f = 1, (G3GV) 10 = 1. 

• Gt = (p 3 ,p 2 ,p 4 ): (GxGV) 15 = 1, (G 2 Gt) 5 = 1, (G3GV) 10 = 1. 

• Gt = (p 2 ,l,p 4 ): (G 1 G e f = 1, (G 2 G,) 15 = 1, (G 3 G,) 10 = 1. 

• Gt = (p 3 , l,p): (G 1 G e ) 5 = 1, (G 2 G e ) 15 = 1, (G3G7) 10 = 1. 

• G £ = (p,p 3 ,p 4 ): (Gift) 5 = 1, (G 2 G £ ) 15 = 1, (G 3 G £ ) 10 = 1. 

• G £ = (p 4 ,p,p 2 ): {G l G l f = 1, (G 2 G £ ) 15 = 1, (G 3 G £ ) 10 = 1. 
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Here p = e 27rl//5 , and we write Gt = diag(a, b, c) by (a, b, c) for notational simplicity. There 
are two sets of solutions: one with m, = (5, 3, 2) (or (3, 5, 2)), and one with m ; = (15, 5, 10) 
(or (5, 15, 10)). Within each set, the two permutations of the rrii are connected because the 
corresponding GYs are related. For example, for the first set of solutions the two options for 
Ge, Ge 1 = diag(l,p, p 4 ) and Ge 2 = diag(l, p 2 , p 3 ), satisfy Gn x = Gf 2 and Gi 2 = G\ . Hence, 
they must generate the same group structure. Similar relations exist for the second set of 
solutions; in addition, the GYs can be obtained by multiplying either GY or Gi 2 by p, p 2 , p 3 
or p 4 and using the freedom to exchange the resulting the (22) and (33) elements of Gi. 

Given that we begin with relations of the form G\ = (GiGe) mi = 1, we seek embeddings in 
the polyhedral group H(m,n), for which the abstract generators s and t satisfy the defining 
relation s 2 = t n = (st) m = 1. Note that H(m,n) = H(n,m) through a redefinition of the 
abstract generators s and t. The number of finite polyhedral groups is strongly constrained, 
as H(m,n) is finite if and only if 2(n + m) > rim [42J. These finite groups are H(n,2), 
#(3,3), #(3,4), and #(3,5), where #(2,2) = # 4 , #(n,2) = D n (n > 3), #(3,3) = A 4 , 
#(3,4) = and #(3,5) = A 5 . 

Here we are considering n = 5, and hence the finite group options are D n (n > 3) and 
A 5 . For the cases with GY and G^ 2 , m ; = 3 and rrii = 2, which corresponds to subgroups of 
A 5 and D 5 , respectively. For rrii — 3 (i.e., the group generated by {Gi, G^} and the group 
generated by {G2, Ge 2 }), the generators obey G\ = 1, G 2 = 1, and (GiGe) 3 = 1. Hence, 
according to Lagrange's theorem of group theory [43], the group order should be divisible by 
30 (the least common multiple of 2, 3, and 5). However, the maximal order of the nontrivial 
subgroups of A 5 is 12, and thus for = 3, the corresponding group is precisely A 5 . By 
similar reasoning, for rrii = 2 the group is D 5 . Since D 5 is a subgroup of A 5 , we conclude that 
the minimal family symmetry group for this class of solutions is indeed A$. It is important to 
note that the groups generated for rrii = 5 (i.e., by {G 2 , GeA or by {G\, Gi 2 }) are the same 
(A§) as the ones generated for the rrii = 3 case given above. For the remaining solutions, 
the fact that the GYs are related to Ge 1 or Ge 2 indicates that these cases generate identical 
groups. The groups generated are either a 300-element group consisting of the elements of 
A5, pA 5 , p 2 A 5 , p 3 A 5 , and p 4 ^4s, or a 50-element group consisting of the elements of D$, pD 5 , 
p 2 D 5 , p 3 D 5 , and p 4 D 5 . Since A 5 is a subgroup of the larger 300-element group, the minimal 
family symmetry group is precisely A$ for all cases. As stated, it is minimal in the sense 
that any other viable horizontal symmetry group should contain A 5 as a subgroup to ensure 
that the full #4 symmetry of the GR1 neutrino mass matrix is preserved. 

We must also show that this is a faithful representation of the group. To this end, we begin 
by expressing the elements of A$ in terms of the generators S and T, which satisfy the A5 
relations S 2 = T 5 = (ST) 3 = 1; the results are shown for completeness in Section 3. Taking 
the specific choice that S = Gi and T = diag(l, p,p 4 ) (or equivalently T = diag(l, p 4 , p))Jj 
it is straightforward to show that all of the representation matrices of the A 5 group elements 
are distinct. This ensures that we indeed have a faithful representation of the group. With 
this choice of S and T, G 2 and G3 are given by 

G 2 = T 3 ST 2 ST 3 , G 3 = T 3 ST 2 ST 3 S. (14) 

Note that we could also choose a representation in which S = G 2 and T = diag(l, p 2 , p 3 ) (or 
4 Note that a similar choice of S and T was taken in |29j . 
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T = diag(l, p 3 , p 2 )), for which G\ and G 3 are then given by Eq. (fT4l) with G 2 — > G±. 

These results were obtained for the choice that n — 5 {i.e., G\ — 1). One might ask the 
question as to whether higher values of n can result in a smaller minimal family symmetry 
group than A%. We have explicitly checked that for a range of n from 6 to 100, the resulting 
family symmetry group is finite if and only if n = 5k, where k is a positive integer. Hence, A§ 
is still the minimal group in these cases. We did not consider the case of n > 100 explicitly, 
but we can argue that the result should hold in the following way. The order j of Gi should 
be a factor of n. However, for j < 100, the horizontal group is finite if and only if j is divisible 
by 5, while if i > 100, the order of the group should also exceed 100. Hence, while it is not 
a mathematical proof, this indicates that the minimal family symmetry group is indeed A5. 
To obtain a more rigorous proof, one can follow the procedure of [44] of enumerating and 
rejecting other relevant finite groups. This topic is left for future work. 

For the case of GR2 mixing, we applied an identical procedure as outlined above based 
on the symmetry generators G\ given in Eq. (1121) . in which we begin by enumerating possible 
choices for n and examining whether such cases can result in a finite family symmetry group 
with one or more of the G'^s. In a systematic scan of the possible values of n up to n = 100, 
we were unable to find solutions for the integers rrti such that (G'^Gg) 171 * = 1, and hence the 
symmetry groups in these cases are infinite. It is worth noting, however, that just as in the 
case of GR1 mixing, one cannot rule out the possibility of a discrete group with a very large 
order. Once again, we defer a more detailed study of this topic for future investigation. Note 
that it can be shown that the dihedral group D\q used in [25] only contains Gg and G' 3 and 
not G[ 2 , and hence it is not minimal according to the criteria specified here. 

For this reason, we now concentrate our attention on the case of GR1 mixing models 
based on the minimal family symmetry group of A 5 for the remainder of this work. After a 
brief summary of group theoretic aspects of A 5 given in the following section of our paper, 
we turn in Section 4 to the question of flavor model building based on A 5 that is predicated 
on our symmetry analysis. The reader that wishes to turn directly to the issue of flavor 
model building should therefore skip Section 3 and proceed to Section 4. 



In this section, we enumerate basic group theoretic properties of A5. A§ is the alternative 
group of order five; i.e., it is the group of even permutations of five objects. From the 
geometrical point of view, A 5 is the subgroup of the three-dimensional rotation group that 
preserves the symmetries of the icosahedron. The order of A 5 is 60; due to this relatively 
large value, the structure of A 5 is more complex than that of many other family symmetry 
groups such as A 4 and S4 that are considered in the literature. The A 5 group elements can 
be generated by two generators S and T that obey the relations: 



3 The discrete group A 5 



S 2 = T 5 = (ST) 3 = 1. 



(15) 



Without loss of generality, we can take (42] 



S=(12)(45), T= (13542), 



(16) 
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in which the cycle (13542) denotes the permutation (1, 3, 5, 4, 2) -> (3, 5, 4, 2, 1), and (12) (45) 
denotes (1, 2, 3, 4, 5) — > (2, 1, 3, 5, 4). The 60 elements of A 5 belong to five conjugacy classes 
d (i — 1, . . . 5). In terms of Schoenflies notation, in which C% labels rotations by 2nk/n and 
the number in front identifies the number of elements in each conjugacy class (and for k — 1, 
the superscript is not written), C\ — 1, C% — 15C*2, C3 = 2OC3, C4 = I2C5, and C5 = 12C 2 . 
The group elements of each conjugacy class can be expressed in terms of S and T as follows: 

Ci : 1, 

C 2 : ST 2 ST 3 S = (12) (34), TST 4 = (13)(24), T A (ST 2 ) 2 = (14)(23), 
T 2 ST 3 = (12) (35), T 2 ST 2 ST 3 S = (13)(25), ST^T = (15)(23), 
5 = (12)(45), T 3 ST 2 ST 3 = (14) (25), T 3 ST 2 ST 3 S = (15)(24), 
T 3 ST 2 = (13)(45), T 4 ST 2 ST 3 S = (14)(35), TST 2 S = (15)(34), 
ST 3 ST 2 S = (23) (45), T 4 ST = (24) (35), (T 2 S) 2 T 4 = (25)(34), 

C 3 : ST = (134), TS = (235), ST 4 = (253), T A S = (143), TST 3 = (145), 
T 2 ST 2 = (243), T 2 ST 4 = (154), T 3 ST = (125), T 3 ST 3 = (234), 
T 4 ST 2 = (152), TST 3 S = (124), T 2 ST 3 S = (354), T 3 ST 2 S = (123), 
ST 2 ST 3 = (345), ST 3 ST = (245), ST 3 ST 2 = (132), (T 2 S) 2 T 2 = (135) 
T 2 (T 2 S) 2 = (254), (ST 2 ) 2 S = (153), (ST 2 ) 2 T 2 = (142), 

C 4 : T = (13542), T 4 = (12453), ST 2 = (14235), T 2 S = (13425), 
ST 3 = (15243), T 3 S = (15324), STS = (12345), TST = (15432), 
TST 2 = (12534), T 2 ST = (14523), T 3 ST 4 = (14352), T 4 ST 3 = (13254), 

C 5 :T 2 = (15234), T 3 = (14325), ST 2 S = (13524), ST 3 S = (14253), 

(ST 2 ) 2 = (12543), (T 2 S) 2 = (14532), (ST 3 ) 2 = (12354), (T 3 S) 2 = (13452), 
(T 2 S) 2 T 3 = (15423), T 3 (ST 2 ) 2 = (15342), T 3 ST 2 ST 4 = (13245), 
T A ST 2 ST 3 = (12435). 

A 5 has 59 subgroups of order 1, 2, 3, 4, 5, 6, 10, 12 or 60, as follows: 

1. The trivial group that consists solely of the identity element. 

2. The two-element subgroup generated by a double transposition: 

Z 2 (1) = {l,ST 2 ST 3 S},Z? ] = {1,TST 4 }, zf^ = {1, T 4 (ST 2 ) 2 }, zf = {1,T 2 ST 3 }, 
Zf = {1, T 2 ST 2 ST 3 S}, Zf = {1, ST 2 ST}, Zp = {1, S}, zf = {1, T 3 ST 2 ST 3 }, 
Z (9) = {i^t 3 ST 2 ST 3 S},z£ 0) = {l,T 3 ST 2 },Zp } = {1, T 4 ST 2 ST 3 S}, 
Z ( 2 12) = {1,TST 2 S},Z { 2 13) = {1,ST 3 ST 2 S},Z^ 14) = {1, T 4 ST}, Z^ 5) = {1, (T 2 S) 2 T 4 }. 

3. The three-element subgroup generated by a 3-cycle: 

Zf ] = {1,T 3 ST 2 S,ST 3 ST 2 },Z^ 2) = {1,TST 3 S,(ST 2 ) 2 T 2 }, zf = {1,T 3 ST,T 4 ST 2 }, 
Zf = {1,ST,T 4 S},^ 5) = {1,(T 2 S) 2 T 2 ,(ST 2 ) 2 S},^ 6) = {1, TST 3 , T 2 ST 4 }, 
Z^ 7) = {1,T 2 ST 2 ,T 3 ST 3 },ZP = {1,TS, ST A },zf ] = {1, ST 3 ST,T 2 (T 2 S) 2 }, 
Zf 0) = {1,ST 2 ST 3 ,T 2 ST 3 S}. 
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4. The four-element subgroup comprising of all double transitions on four of five objects, 
which is isomorphic to K±: 



K 
K 
K 



P = {l,ST 2 ST 3 S,TST A ,T\ST 2 ) 2 },Kf ] 
f = {l,S,T 3 ST 2 ST 3 ,T 3 ST 2 ST 3 S},K^ 
(5) - {1, ST 3 ST 2 S, T 4 ST, (T 2 S) 2 T 4 }. 



{1, T 2 ST 3 , T 2 ST 2 ST 3 S, ST 2 ST}, 
{1, T 3 ST 2 , T 4 ST 2 ST 3 S, TST 2 S}, 



5. The five-element subgroup generated by a 5-cycle, which is isomorphic to Z^. 



r(l) 

'5 

,(3) 
'5 

,(5) 



{1, STS, ST 2 S, ST 3 S, TST}, zf ] = {1, ST 3 , T 2 S, {ST 3 ) 2 , (T 2 S) 2 } 



3\2 



{1, T 2 ST, T 4 ST 3 , T 3 (ST 2 ) 2 , T 4 ST 2 ST 3 }, zf ] = {1,T, T 2 , T 3 , T 4 }, 



{1, TST 2 , T 3 ST\ (T 2 S) 2 T 3 , T 3 ST 2 ST 4 }, zf = {1, ST 2 , T 3 S, (ST 2 ) 2 , (T 3 S) 2 }. 



6. The subgroup of order six in which three objects obey symmetric permutations, which 
is isomorphic to S3: 



(10) 



{1,S, T 3 ST 2 , T 3 ST 2 S, ST 3 ST 2 , ST 3 ST 2 S}, 
{1,T 2 ST 3 , T 4 ST, TST 3 S, {ST 2 ) 2 T 2 , T*ST 2 ST 3 S}, 
{1,T 3 ST, T 4 ST 2 , TST 2 S, (T 2 S) 2 T A , ST 2 ST 3 S}, 
{1, ST, T 4 S, {T 2 S) 2 T\ T 3 ST 2 ST 3 , T 2 ST 2 ST 3 S}, 
{1, TST 4 , T*ST, (T 2 S) 2 T 2 , {ST 2 ) 2 S, T 3 ST 2 ST 3 S}, 
{1, TST 3 , T 2 ST\ ST 2 ST, T 4 (ST 2 ) 2 , ST 3 ST 2 S}, 
{1,T 2 ST 2 , T 3 ST 3 , ST 2 ST, TST 2 S, T 3 ST 2 ST 3 S}, 
{1, TS, ST 4 , T 4 (ST 2 ) 2 , T 3 ST 2 ST 3 , T 4 ST 2 ST 3 S}, 
{1, TST 4 , T 3 ST 2 , ST 3 ST, T 2 (T 2 S) 2 , T 2 ST 2 ST 3 S}, 
--{1,S, T 2 ST 3 , ST 2 ST 3 , T 2 ST 3 S, ST 2 ST 3 S}. 



7. The order ten subgroup generated by a 5-cycle and a double transition which conjugates 
it to its inverse, which is isomorphic to D 5 : 



D 
D 
D 
D 
D 
D 



{1, STS, ST 2 S, ST 3 S, TST, T 2 ST 3 , T 3 ST 2 , T 4 (ST 2 ) 2 , (T 2 S) 2 T 4 , T 3 ST 2 ST 3 S}, 
{1, ST 3 , T 2 S, (ST 3 ) 2 , (T 2 S) 2 , T 3 ST 2 , T 4 ST, ST 2 ST, ST 2 ST 3 S, T 3 ST 2 ST 3 }, 
{1, S, TST 4 , T 2 ST, T 4 ST 3 , ST 2 ST, T 3 (ST 2 ) 2 , (T 2 S) 2 T 4 , T 4 ST 2 ST 3 , T 4 ST 2 ST 3 S}, 
{1, T, T 2 , T 3 , T 4 , ST 2 ST 3 S, ST 3 ST 2 S, T 2 ST 2 ST 3 S, T 3 ST 2 ST 3 S, T A ST 2 ST 3 S}, 
{1, S, TST 2 , T 3 ST 4 , T 4 ST, TST 2 S, (T 2 S) 2 T 3 , T 4 (ST 2 ) 2 , T 3 ST 2 ST\ T 2 ST 2 ST 3 S}, 
{1, ST 2 , T 3 S, (ST 2 ) 2 , (T 3 S) 2 , TST 4 , T 2 ST 3 , TST 2 S, ST 3 ST 2 S, T 3 ST 2 ST 3 }. 



8. The twelve-element subgroup of permutations of any four of the five objects, which is 
isomorphic to A A : 
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Conjugacy Classes 


Ci 


c 2 


c 3 


c 4 


c 5 


1 


1 


1 


1 


1 


1 


3 


3 


-1 









3' 


3 


-1 









4 


4 





1 


-1 


-1 


5 


5 


1 


-1 









Table 2: The character table of the A$ group, in which the conjugacy classes written in Schocnflics notation 
are d = 1, C 2 = 15C 2 , C 3 = 20C 3 , C 4 = 12C 5 , and C 5 = 12Cf . 



4 1} = {1, ST, T A S, TST 4 , T 2 ST 2 , T 3 ST 3 , ST 3 ST 2 , TST 3 S, T 3 ST 2 S, (ST 2 ) 2 T 2 , T\ST 2 ) 2 , 
ST 2 ST 3 S}, 

Af = {1, ST\ TS, T 2 ST 3 , T 3 ST, T 4 ST 2 , ST 2 ST, ST 3 ST 2 , T 3 ST 2 S, {T 2 S) 2 T 2 , (ST 2 ) 2 S, 
T 2 ST 2 ST 3 S}, 

4 3) = {1, S, TST 3 , T 2 ST\ T 3 ST, T 4 ST 2 , ST 3 ST, TST 3 S, (ST 2 ) 2 T 2 , T 2 {T 2 S) 2 , T 3 ST 2 ST 3 , 

T 3 ST 2 ST 3 S}, 

4 4) = {1, ST, T 4 S, TST 3 , T 2 ST\ T 3 ST 2 , ST 2 ST 3 , TST 2 S, T 2 ST 3 S, {ST 2 ) 2 S, {T 2 S) 2 T 2 , 

T*ST 2 ST 3 S}, 

Af = {1, ST\ TS, T 2 ST 2 , T 3 ST 3 , T A ST, ST 2 ST 3 , ST 3 ST, T 2 ST 3 S, T 2 {T 2 S) 2 , {T 2 S) 2 T\ 
ST 3 ST 2 S}. 



9. The full group. 

For completeness, we now summarize a number of basic group theoretic facts about A 5 that 
are needed for flavor model building, which have previously been presented in the literature. 
We begin by recalling that A 5 has five irreducible representations, 1, 3, 3', 4 and 5. The 
A 5 character table is presented in Table [21 and the products of irreducible representations, 
which are easily determined from the character table, are listed below: 

1®R = R®1 = R, 3®3 = 1©3©5, 3'®3' = 1©3'©5, 3x3' = 4© 5, 
3©4 = 3'©4©5, 3'©4 = 3©4©5, 3®5 = 3©3 / ©4©5, 
3' ©5 = 3© 3' © 4 © 5, 4©4 = 1©3©3'©4©5, 4©5 = 3©3 / ©4©5i©5 2 , 
5©5 = 1©3©3' ©4ffi4ffi5 1 ©5 2 . (17) 

Here R indicates any A 5 irreducible representation and 5x,2 stand for the two 5 representa- 
tions that appear in the Kronecker products. In the Appendix, we list the explicit represen- 
tation matrices of the generators S and T for each of the irreducible representations, and 
provide the corresponding Clebsch-Gordan coefficients. 
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4 A 5 Flavor Models with Triplet Embeddings 



Guided by the above symmetry analysis, we now turn to the issue of lepton flavor model- 
building based on A5. Lepton flavor models have previously been constructed within this 
framework [271129] . In this work, we specifically seek natural models of GR1 mixing in 
which the SM leptons are embedded in the triplet representations (as opposed to singlet 
representations, as done in [29]), and in which issues of vacuum alignment can be achieved 
naturally. In this section, we will present two A§ flavor models that result in GR1 mixing: 
one in which the family symmetry is broken to Z5 and K4 in the charged lepton and neutrino 
sectors, respectively, and one in which the Z 5 symmetry is itself broken in the charged lepton 
sector. The details of each scenario will be described below. However, we begin by reviewing 
basic aspects of flavor model building within the framework of A 5 . 

The basic paradigm is that the full family symmetry under which the matter fields trans- 
form is realized at some high energy scale A, and is then spontaneously broken to subgroups 
by flavon fields, which are scalar multiplets with nonvanishing vacuum expectation values 
(VEVs). For the leptons, the subgroup that is preserved in the charged lepton mass matrix 
should be different than the surviving subgroup in the neutrino sector. This misalignment 
between the symmetry breaking patterns results in mass-independent textures (such as GR1 
or tri-bimaximal mixing) at leading order (LO). Hence, group theoretical considerations of 
the family symmetry breaking chains are critical for guiding model-building efforts. 

As we have demonstrated in section [21 the GR1 mixing pattern is naturally produced 
in A 5 models when the A 5 symmetry is broken to the Z 5 subgroup generated by T in the 
charged lepton sector, and broken down to the K4 subgroup comprised of S, T 3 ST 2 ST 3 , 
T 3 ST 2 ST 3 S and the identity in the neutrino sector. Since the lepton masses are generated 
after the spontaneous breaking of the flavor symmetry, the flavon VEVs must be invariant 
under the corresponding subgroups, which strongly constrains the vacuum alignment in the 
flavon field sector. More precisely, for the flavon field ipx associated with the charged leptons 
at LO, we see from the explicit representation matrices shown in the Appendix that the 
invariance under the action of T means we should restrict ourselves to the following set of 
flavon field VEVs: 



(<Pt) oc 



(1,0,0), <^ T ~3 

(1,0,0), <^ T ~3' 

(0,0,0,0), <^ T ~4 

(1,0,0,0,0), ifr^h. 



Similar reasoning for the flavon field <ps associated with the K4 subgroup of the neutrino 
sector, its vacuum configuration is determined to be 

r (1,1,1, ik ^ 5 ~4 

in which v, v 2 and V3 are arbitrary parameters. It is straightforward to see that if ips is a 
A§ triplet 3 or 3', there is no vacuum alignment that maintains the K4 symmetry, although 
such representations can break A 5 to the Z 2 subgroup generated by S and T 3 ST 2 ST 3 or 
T 3 ST 2 ST 3 S separately. 
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Fields 


L 


E c 


N c 


h u ,d 




X 





A 


<P 


X° 


0° 


A u 


^° 




3 


3 


3' 


1 


1 


3 


5 


4 


5 


3 


5 


4 


3' 


z 5 


P 3 


1 


p 2 


1 


P 2 


P 2 


p 2 


1 


P 


P 


P 


1 


P 4 



Table 3: The transformation properties of the matter fields, the electroweak Higgs doublets, the flavon fields 
and the driving fields under the flavor symmetry A§ x Z5, where p = e 27 ™/ 5 , L = (^1,^2,^3) are the lepton 
doublet, E c = (e c , r c , /j, c ) is the charged leptons, and N c = (vf, i/f, v§) are the right-handed neutrinos. 



4.1 Model 1 

Armed with this knowledge, we now present the first model (Model 1), in which the full 
flavor symmetry is A 5 xZ 5 . The particle content is shown in Table [3J In keeping with our goal 
of triplet embeddings, the lepton doublets L = (^1,^2,^3) and right-handed (RH) charged 
leptons E° = (e c , r c , /i c ) are both assigned to A 5 triplet (3) representations. We note that r c 
and /i c are taken to be the second and third components respectively so that the LO charged 
lepton mass matrix is diagonal, as will be seen below. The neutrino masses are generated 
via the Type I seesaw mechanism, in which the right-handed neutrinos N c = v 1 ^, ^3) 
transform as a 3'. The flavon sector consists of the following set of fields: a singlet £, triplet 
X, and a 5-plet (f>, which couple to the charged leptons, and the 4-plet A and additional 5-plet 
tp, which couple to the neutrino sector. We use the now-standard supersymmetric driving 
field method to arrange the vacuum alignment, in which a continuous U(1)r symmetry 
related to R— parity is introduced under which the matter fields have a +1 R— charge while 
the electroweak Higgs and flavon fields are uncharged. Upon the introduction of so-called 
driving fields that carry a U(1)r charge of +2, we see that the flavon sector superpotential 
self-interactions must consist of terms in which the driving fields appear linearly. Hence, the 
minimization of the flavon potential can be achieved simply by ensuring that the F-terms of 
the driving fields vanish at the minimum. In our model, we require four driving fields: two 
triplets, one 3 (%°) and one 3' (V ;0 ), one 4 (A ), and one 5 (0°). With this particle content, 
the superpotential w takes the following form at leading order: 

w = wi + w y + w d , (20) 

in which 

Wi = V -^{E c LUh d + y^((E c L) 3X hh d + ^((M) 5 0)i^ (21) 

w v = ^((7V c L) 4 A) 1 /i u + ^x((iV c iV c )5^)i (22) 

w d = A(X°XH + / 2 (x°(x0) 3 )i + Si(0V)i£ + 92(<t>°(xxh)i + S 3 (0°(x0) 5 )i (23) 
+ #4(0 o (#) 5l )i + & 5 (0°(#)5 2 )i + M A (A°A) 1 + /n(A°(AA) 4 )i + /i 2 (^°(A^) 3 0i, 

in which the subscripts 1, 3, 3', 4 and 5 denote A 5 contractions. In the above, we have 
neglected to write higher dimensional operators, which are suppressed by additional powers 
of the high (cutoff) scale A. We note that the combinations (x°(xx)s)i an d (x° (fid) 3)1 are 
omitted here since the contractions {xx)s an d (00)3 vanish. 
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The F-terms of the driving fields take the form 

^ = hxit + / 2 (-2xi0i + Vz + Vs X30 2 ) = o 

= /lX3^+ /2(V3Xl^5 - ^^204 + X30l) = 
q-q = /lX2^+ 72(^X102 + X201 - ^^303) = 

-qJP = 9l<h.t + ~ X2X3) + V / 3^ 3 (X205 - X302) + 2^4(0? + 0205 ~ 20 4 4 ) 

+2g 5 (<pl - 20 2 05 + 3 4 ) = 

-Qjp = gifeU - 92XlX3 + gsiXlfe + ^2 X204 + X30l) + 2# 4 (0105 + V& 0204) 

+^(-40105 + ^0^) =0 

-Qjjt = £i0 4 £ + VQg2xl + ^s(2xi04 + v / 2x30s) + ^ 4 (-40i0 4 + VQ(pl) 

+ 2^5(0104 + ^0203) =0 

^4 = 01<fof + ^£2X2 + fi'3(-2Xl03 - ^ ^202) + fl'4(-40103 + \f§ 4>l) 
O04 

+ 2^(0x03 + ^0405) =0 
^To = ^102^ - 2^^2X1X2 + g3(-Xl<l>2 ~ VSX2(f>l - ^2 X303) + 2^(0102 + VQ 3 5 ) 

+^(-40102 + ^0^) = 0. (24) 



A solution to these relations is: 

(0 = ^> (x> = K,o,o), (0) = (^,0,0,0,0), (25) 
in which the VEVs obey the relations 

/1 

u = 7Tr~ v £i QiVjOt + 25f 2 f J + 2(#4 + £ 5 )^ = 0, (26) 

^/2 

with undetermined. The vacuum expectation values of A and tp, which give rise to GR1 
mixing in the neutrino sector, are determined by the F-terms of the associated driving fields 
as follows: 

9W d _ A , . /OA A , A2^ 



<9A° 

dw d 
dA° 2 
dw d 
dA° 3 
dw d 
dA° 4 



M A A 4 + /n(2A 1 A 3 + A^) = 
M A A 3 + /i 1 (2A 1 A 2 + A*) = 
M A A 2 + fci(A? + 2A3A4) = 
M A Ai + /ii(2A 2 A 4 + Aij) = 

13 



dw, 
dwd 
dm 

which admit the solution 



V2h 2 (A 1 ip 5 + 2A 2V 9 4 - 2A 3 (^ 3 - A 4 ^ 2 ) = 
h 2 (-2A l( p 3 + A 2 <^ 2 + VQ A- m - 3A 4 ^ 5 ) = 
h 2 (3A l (p 2 - VQA 2 (pi - A 3 ip 5 + 2A 4 (p A ) = 0, 



(A) = (1, 1, 1, l)v A , (y?) = (\l -(v 2 + u 3 ), v 2 , v 3 , v 3 , v 2 ) 



(27) 



(28) 



where v& = —M&/(3hi) and v 2y3 are undetermined. The auxiliary symmetry Z5 ensures 
that the flavons £, x an d 4> which appear in u>£ don't contribute to w v and vice versa at 
LO. With this set of flavon VEVs, A 5 is broken to Z 5 and K± in the charged lepton and 
neutrino sectors, respectively. After flavor and electroweak symmetry breaking, we results 
in the following diagonal charged lepton mass matrix: 












U0 












(29) 



in which Vd = (hd) is the VEV of the electroweak Higgs field hd- The three independent 
combinations y^v^Vd/A, ye 2 v x Vd/A and yg 3 v^Vd/A can be expressed in terms of the charged 
lepton masses as follows: 



Vix-^d 



1 

t 

1 



2 

v 4> 1 to \ 
ye 3 -£V d = -{2m e - m M - m T ). 



(30) 



We see that y^v^ : ye 2 v x : y^v^ ^2:3: —1. Hence, to reproduce the observed mass 
hierarchies of the charged leptons, additional fine tuning of the parameters ye t is required. 

Turning now to the neutrino sector, the Dirac neutrino mass matrix as obtained from 
the first term of Eq. (1221) is given by 



m D = y u 




(31) 



in which v u = (h u ). The second term in Eq. 022|) results in the Majorana mass matrix 

( 2^\{v 2 + v z ) ~V3v 3 -y/3v 3 \ 

-V3v 3 \/Qv 2 -\\(v 2 + v 3 ) 



\ 



V3 



''3 



■S («2 + V3) 



''2 



(32) 



/ 
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The light neutrino mass matrix m„ = —nv^m^niD is exactly diagonalized by the GR1 
mixing matrix Ugri via ^qri^^gri = diag(mi, m 2 , m 3 ), in which £/gri is given by Eq. (jlj) 
with cot 812 = 4> g . The mass eigenvalues are 

4:V6ylv 2 A vl AVGylv 2 A vl 1^y 2 v v\vl 

(a - 6) A A (a + 6) A A aA A v 7 

in which the parameters a and 6 are defined as 

a = x (Av 2 + v 3 ), b = 3V5xv 3 . (34) 

The neutrino mass spectrum is strongly constrained in the present model, as it is determined 
by three independent real parameters: the absolute values of a and b and the relative phase 
between them. In particular, note that the following sum rule is obeyed: 

111 
— + — = —. 35 

rrii m 2 ^3 

This sum rule has also been obtained in a number of other models [UJ . The mass spectrum 
can either be a normal hierarchy (NH) or inverted hierarchy (IH). Taking into account the 
measured solar and atmospheric mass-squared differences, we obtain the following limits for 
the lightest neutrino mass: 

\nii\ > 0.011 eV, NH 

\m 3 \ > 0.028 eV, IH. (36) 

We note that an analytical formula for the lower limits of neutrino masses that obey this 
sum rule has been determined in (UJ. In addition, this model has definite predictions for the 
Majorana phases ai 3 and 023 (using the standard parameterization [45] of the MNSP mixing 
matrix). These phases are determined by the phase differences between the complex masses 
of Eq. ( 1331) . After some algebraic manipulation, the following relations can be obtained: 

costti 3 = 1, 1 -[|mi| 2 |m 2 | 2 - (|mi| 2 - |m 2 | 2 )|m 3 | 2 l 

2|mi||m2||J7i3| 

COSa 23 = 771 PJT n f[l m l| 2 | m 2| 2 + (|«1l| 2 - |«l2| 2 )|"l3| 2 ]- (37) 

2|mi| z |m2| |77i3| 

Finally, we comment that in this model, the neutrinoless double beta decay parameter m ee 
can also be expressed in terms of the light neutrino masses and the golden ratio (f) g as follows: 

1 2 1 |2 , i I™ 1 2 1 |2 1 I™ 1 2 1 I21 1 / 2 



\m P 



x/5 



m 3 \ 



[\nii\ 2 \ni2\ 2 + (f)g\nii\ 2 \m 3 \ 2 — —\ni2\ 2 \m 3 \ 2 ^ . (38) 



The lower bound of \m ee \ is found to be about 4 meV and 52 meV for the NH and IH cases, 
respectively. While the GR1 mixing pattern emerges naturally within this model, we must 
tune the Yukawa couplings to fit the charged lepton masses. Further fine tuning is re- 
quired if subleading corrections are included. We note that this kind of fine tuning does not 
contradict the logic of Section 2, which only focuses on flavor mixing. In what follows, we 
aim to improve on this scenario by proposing another model that overcomes this defect 
by breaking the residual Z 5 symmetry in the charged lepton sector. 
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4.2 Model 2 



4.2.1 LO predictions 

To improve on the previous model, note that there are four possible assignments for the 
lepton doublets and the charged leptons: (i) both can be 3 or 3' representations, (ii) one 
is a 3 and the other is a 3', (iii) one is a 3 or 3' while the other is an A 5 singlet, and (iv) 
both are A$ singlets. We set aside case (iv), in which does not play a fundamental role in 
the lepton sector. In case (ii), which was considered in [27], it is in general difficult to avoid 
nonvanishing off-diagonal elements in the charged lepton mass matrix, which will result in 
deviations from GR1 mixing. Case (iii) was considered by [29], in which an additional Z n 
was introduced to ensure a diagonal charged lepton mass matrix; in this model, the necessary 
vacuum alignment is not realized consistently within the driving fields method. 



Fields 


L 


E c 


N c 


h u ,d 


X 





c 




e 


p° 








3 


3 


3 


1 


5 


5 


1 


5 


i 


3' 


5 


5 




1 


UJ 2 


1 


1 


UJ 2 


UJ 


1 


1 


i 


1 


UJ 


1 


z 3 


UJ 2 


UJ 


1 


1 


1 


1 


UJ 


UJ 


i 


1 


1 


UJ 



Table 4: The transformation properties of the matter fields, the Higgs doublets, the flavon fields and the 
driving fields under the flavor symmetry A5 x Z% x Z3, in which uj = e 2lri / 3 , L = (^1,^2,^3) are the lepton 
doublet, E c = (e c ,^ c ,r c ) are the charged leptons, and N c = {y\,v%,v%) are the right-handed neutrinos. 



Given that we wish to focus on triplet embeddings that result in GR1 mixing without 
excessive fine-tuning, we thus focus on case (i), and assign the SM leptons to 3 representations 
(3' representations give essentially the same results). The family symmetry of our second 
model is A^ x Z% x Z3, and the particle content is given in Table HI We see that in this case 
the flavon fields are only in 1 (£) and 5 (x, 4>, V 9 ) representations. Four driving fields are 
needed: one singlet one 3' (p°), and two 5's (x° and if ). As before, the superpotential 
of the model takes the form 

w = we + w v + w d . (39) 
Let us focus first on the driving superpotential Wd, which at LO in this model is given by 
w d = <?i£°(x0)i + 92{ P \xMx + M x ( X °x)i + S3(X°(00) 5l )l + S4(X°(00)5 2 )l 

+ MAOiC + ^(^W) 5l )i + /* 3 (</(w)5 2 )i- (40) 

The charge assignments reported in Table H] indicate that the minimization equations of the 
driving potential can be separated into two decoupled sets: one for the neutrino sector and 
one for the charged lepton sector. In the charged lepton sector, the minimum is determined 
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by the following conditions: 

dw d 



S'l(Xl01 + X205 + X304 + X403 + Xbfa) = 
#2(2X205 - X304 + X403 ~ 2X502) = 
5'2(-V / 3Xl04 + V / 2~X203 - V / 2~X302 + V / 3~X40l) = 
#2(^X103 - V / 3X301 + V / 2~X405 - ^^504) = 
M xX l + 2# 3 (0? + 0205 - 20 3 4 ) + 2# 4 (0? - 20 2 5 + 0304 
M X X5 + 2^(0105 + ^0204) + # 4 (-40!0 5 + VQ <j>\) = 
M x Xi + ^("40104 + V§4>1) + 2^(0104 + VQ 0203) = 



dwa 
dpi 
dm 
dpi 
dwd 
dpi 
dwd 

d X °i 
dwa 

dxl 
dwd 
dxl 

= M xX 3 + ^("40103 + V64>1) + 2^(0103 + ^0405) = 

dw, 



dxl 



M xX 2 + 2(?3(0102 + ^0305) + ^(-40102 + V^^) = 0, (41) 
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which result in the solutionis 

( X ) = (0, 0, 0, 0, v x ), (0) = (0, 0, v h 0, 0), (42) 

in which 



-VQ g4V^/M x , undetermined. (43) 

Under the action of T, we have T(x) = p 4 (x) an d T(</>) = p 2 (0). As a result, the Z 5 
symmetry is broken completely: no residual symmetry is preserved in the charged lepton 
sector. For the neutrino sector, the minimization equations take the following form: 



dwd 
<M 
dwd 
d<p° 
dwd 

dw d 

art 

dwd 
d<p° 



hiifiC + 2h 2 {(f\ + (p 2 ip 5 - 2^4) + 2h 3 (ipl - 2(p 2 (p 5 + ^3^4) = 

hxipnC + 2/i 2 (y?iy?5 + Vg^'Pa) + h 3 (-iipiip 5 + v6 ipf) = 

htpiC + h 2 (-4:Lp 1 Lp i + V(i<pl) + 2h 3 ((p 1 (p 4 + V6(p 2 ip 3 ) = 

hx(p 3 ( + h 2 (-4:Lp 1 (p 3 + VQ(p 2 2 ) + 2/i 3 (v9iv?3 + y/E tp^s) = 

h(p 2 ( + 2/i 2 (y?iy?2 + V6ip 3 ip 5 ) + /i 3 (-4^i<£> 2 + v 7 ^ y^) = 0. (44) 



5 There is another solution in which (x) = (0, v x , 0, 0, 0), and (</>) = (0, 0, 0, u^, 0), with w x = —y/6g 4 v?/M x 
and «0 undetermined. However, this solution can be obtained by acting on Eq. (|42[) with the A 5 group element 
T 3 ST 2 ST 3 S , and consequently it is an equivalent solution. 
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The minimum thus takes the form 



(C)=^o (v) = (v ^(v 2 + v 3 ),v 2 ,v 3 ,v 3 ,v 2 ), (45) 



in which 



2\ 



h 2 v 2 (-v 2 + 2w 2 f3 + 4t> 3 ) = h 3 v 3 (-v 3 + 2t> 2 t;3 + 4w 27 
[2h 2 , . [2 h 3r 3v^ , . -1 



In this sector, A 5 is broken to its K± subgroup by the nonvanishing VEVs of ( and tp. Given 
the symmetry of w d , we can generate other minima of the scalar potential by acting on 
Eq. (I42p and Eq. (145]) with group elements. However, these new minima are related to the 
original minimum by field redefinitions, and hence are physically equivalent. As a result, 
there is no loss of generality in selecting Eq. ( )42l) and Eq. ( 1451) as the local minimum. 

We turn now to the Yukawa superpotential for the charged leptons, which in this model 
takes the form 

wt = ^({ETLhfahi + j^(E c L)i(xx)ih d + V -^{{E c L) s {xxh)xh d 

+ V -^{{E c LUxx)^)xh d + ^((E'LUxxhM* + ■■■, (47) 

in which we have neglected to write higher dimensional operators that will be addressed later 
in the paper. After electroweak breaking and flavor symmetry breaking via Eq. ( 1421) . we see 
that the charged lepton mass matrix is diagonal: 



m t = v d . (4 





Clearly, the electron is massless at leading order; it will be shown later that its mass arises 
at next to leading (NLO) order. As the tau and muon masses are suppressed respectively by 
1/A and 1/A 2 , the observed mass hierarchies can be reproduced naturally if v^/A and v x /A 
are of order A 2 , in which A c ~ 0.22 is the Cabibbo angle. We note that the charged lepton 
mass hierarchies can thus be determined by the flavor symmetry itself without invoking 
further Froggatt-Nielsen mechanisms or fine tuning the couplings. As is common in discrete 
flavor symmetry model building, we assume that all VEVs are approximately at the same 
scale. Hence we take v^/A ~ v x /A ~ v^/A ~ v 2 /A ~ v 3 /A ~ in the following. 
In the neutrino sector, the superpotential w v is given as follows: 

w v = ^{N c LUh u + ^{{N c L)w) x K + ^M(V C V C ) X + ... (49) 

where once again we have neglected at this stage to write higher-order contributions. Note 
that we can always choose the mass parameter M to be real and positive by performing a 
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global phase transformation on the right-handed neutrinos. Given the set of flavon VEVs of 
Eq. (I4"5|) . the neutrino Dirac mass matrix is 



/ 



f2 y i^ 2 (>2 + v 3 ) 

-V3y„ 2 v 2 



\ 



-V3y U2 v 2 
VQy^ 2 V3 



y^ 



y^ v c~ \/f y^( v 2 + v3.) 

and right-handed Majorana mass matrix is 



V^yu 2 V2 

yu 2 (v 2 + v 3 ] 
VQyv 2 v3 



\ 



f,(50) 



m M 



MOO 
M 
M 



(51) 



Note that the three right-handed neutrinos are exactly degenerate. Integrating out the heavy 
degrees of freedom, we obtain the light neutrino mass matrix: 



m v = -m T D m h }m D = [/*diag(mi, m 2 , m 3 )£^. 
In the above, the unitary matrix U v is given by 

U v = f/ GR idiag(ze- lQl ,ze- ia2 ,e- ia3 ), 



(52) 



(53) 



in which £/gri is given by Eq. (TjJ for the GR1 mixing case (cot d\ 2 = 4> g ), and the phases 
are given by 



«! = arg(A + B + C), a 2 = arg(A + B - C), a 3 = arg(A - 25), 



with 



5 



2/^ v 2 + 4u 3 



v/6 A 

The light neutrino masses mi, m 2 and m 3 are thus 



C 



30^^ 
A 



(54) 



(55) 



mi = \A + B + C\ 2 ^, 



m 2 = \A + B-C\^, 



m 3 



1 1 M 



(56) 



The neutrino masses in this model thus depend on three complex parameters A, B and C 
as well as the overall factor v\jM. Therefore, mi j2j3 a priori are unrelated and any type of 
mass hierarchy can be accommodated. In particular, there is no general neutrino mass sum 
rule constraint in this case. 



4.2.2 NLO corrections 

In this section, we address the next-to-leading corrections to the results presented above 
for our A 5 x Z 3 x Z 3 model. Let us begin with the driving superpotential. Taking into 
account the NLO operators, 

w d = w° d + 5w d , (57) 
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in which w° d is the leading order contribution presented previously, and Swd is the most 
general A 5 x Z 3 x Z 3 invariant quartic polynomial that is linear in the driving fields. More 
precisely, 5wd can be written as 

1 8 1 19 1 11 

to, = I £ z.xf ^ + I £ r,xf + I £ cWf , (58) 

i=l i=l i=l 

where the 0(1) coefficients a?j, rj and Cj are unconstrained by the family symmetry, and if , 
Xf and If represent a basis of independent quartic invariants, in which the l\ are 



xf = ^°(x(xx)5 1 )i, xf° = e°(x(xx)5 2 )i, xf° = eW^sji, 4° = e° 



5 2 n> 



xf = e°c 3 , xf = e°c(w)i, xf = e°(^(w) 5l )i, 4° = e°(^(w) 52 ) 1 , (59) 

while the Xf take the form 

xf = (p o (x(xx) 3 ) 3 0i, xf = (AX(XX)3') 3 ')1, = (p o (x(xx)4 s )30i, 
Xf = (p°(x(xx) 4 J 3 0i, Xf = (p (x(xx) 5l )3')i, xf = (p°( x (xx)Mx, 
Xf = (p°(0(#)a)3')i, Xf = (p°(0(#)aO8Oi, Xf = (p o (0(#) 4s ) 3 Oi, 

Xf = (p°(H^Kh)l, Xf = (p O (0(00) 5l )3Ol, Xf = (p°(<K#) 52 ) 3 ')l, 

Xfs = (p°(<P<p)a>)iC, Xf = (p°(<p((p(p)s)a>)i, Xf 5 = (p°((p(<p<p)a') 3 ')i, 

Xfe = (p (<p(w)4 s ) 3 ')i, Xf 7 = (p°(^(<p<p) 4A ) 3 ')i, Xfg = (p°((p( W ) 5l ) 3 ')i> 

Xf = (p o (^( W ) 52 ) 3 0i, (60) 

o 

and the Xf are given by 

xf = (xV)i(xx)i, Xf = (x°(0(xx) 3 )5) 1 , if = (X°(0(XX)3')5)1, 

xf = (x o (0(xx)4 S ) 5l )i, xf = (X°(0(XX)4 S ) 52 )1, xf = (X°(0(XX)4J 51 )1, 

Xf = (X°(0(XX) 4 J 52 ) 1 , Xf = (xWxxkk)!, xf = (X°(0(XX) 51 )5 2 )1, 

Xf = (x°(<f>(xx)Mi, Xf = (X°(0(XX) 52 )5 2 )1. (61) 

Note that due to the auxiliary Z 3 x Z 3 symmetry of our model, the subleading corrections 
to the driving superpotential that are proportional to tp° are suppressed by 1/A 2 . Therefore, 
the solutions for the minima of £ and tp given in Eq. (|45]1 is unchanged at this order and it 
is only modified at the next to next to leading order (NNLO). For the remaining fields, the 
new set of vacuum values is obtained as usual by solving the vanishing F-term constraints. 
In so doing, we seek a solution that perturbs the original set of minima to first order in 1/A. 
After a lengthy calculation, we find that there are nontrivial shifts of and x- 

(x) = (8v Xl , 0,5v X3 ,5v X4 ,v x ), {(/>) = (Sv^SvfaVfrSvfaO). (62) 
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Note that the VEVs of the second component of x an d the fifth component of still vanish. 
For the remaining fields, the shifts Sv xi , Sv X3 , 5v X4 , Sv^, 5v^ 2 , and 5v ( f >4 are given as follows: 



5v Xl = [V2d 3 g 2 + 2d 2 (g 4 - 2g- c 

(-Ag 3 + 2g 4 )v x 
5v xs = 7^ 5v 4>i 

V 6 g^Vtf, 



X 



4V3g 2 g 3 Av<p 



Sv X4 



8v 



cj}2 



2d^l 

3#iA^ 

div 3 
3g±Av 



x 



v: 



5v^ 4 = (V2 d 2 g 4 + d 3 g 2 ) 4 ^ , ( 63 ) 

in which the parameters di, d 2 and d 3 are 

di = x 5 + Ax & (2v 2 2 + v 2 v 3 + 2v 2 3 )/(3v 2 ) + (2/3) 3/2 x 7 (v 2 + 4v 3 )(llv| - 2v 2 u 3 - 4w 3 2 )/^ 

+ (2/3) 3 / 2 x 8 (4t; 2 + v 3 )(llv 2 3 - 2 ^3 - 
d 2 = -2v / 3(2r 3 + r 5 ) 

d 3 = -AVQc 4 + 1QV3c 5 -2VQc 8 + VQc 9 . (64) 

Parametrizing the ratio of a generic flavon VEV to the cutoff scale A by the small quantity 
s, we see that 5v xi4 /v x and 8v^ 24 /v^ are both of 0(e). Furthermore, Sv^, which contributes 
to the mass of the electron, is undetermined at this order. To obtain a phenomenologically 
acceptable range of values for the electron mass, we therefore require that Sv^ is O^v^). As 
a result of this choice, 5v X3 is suppressed by e 2 with respect to v x as well. This is clearly an 
assumption in this model, since in the absence of a theory for the higher-order terms, we have 
allowed for the most general higher-order corrections. Indeed, it is a result that ultimately 
is due to our insistence on triplet embeddings of both the left-handed and right-handed SM 
leptons. It is worth noting that this issue can in principle be overcome without resorting to 
triplet embeddings by extending the flavor symmetry to its double cover X', which has 
doublet representations in which the right-handed charged leptons can be embedded (note 
that X' has previously been used in this way for quark flavor model building [30]). 

We turn now to the NLO corrections to the lepton masses and mixings. In the charged 
lepton sector, the following operators can contribute: 

(E c V) x (x(H)^)xh<ilk\ (E c L) 1 ( X ((f>4>) 52 )xh d /A 3 , ((£ c L) 3 (x(00) 3 ) 3 WA 3 , 
((E c L) 3 (x(Whh)ih d /A 3 , ((^ c L) 3 (x(#)4 s ) 3 )iVA 3 , ((M) 3 (x(#) 4 J 3 )iVA 3 , 
((E c LUx(<P<P)Mih d /k\ ((E c L) 3 (x(<P<Ph 2 h)ih d /A 3 , ((E c L) 5X )i(H)ih d /A 3 , 
((E c L) 5 ( X (00) 3 ) 5 ) 1 / id /A 3 , ((E c L) 5 ( X (00) 3 O5)iWA 3 , ((E c L) 5 ( X (00) 4s ) 5l ) 1 VA 3 , 
((M) 5 (x(#)4 s ) 52 )iVA 3 , ((E^Uxm^hJihd/A 3 , ((E c L) 5 ( X (H)4 A h 2 )ih d /A 3 , 
((E^LUxi^hMA 3 , ((E^Uxi^hJ^hJA 3 , {{E c LUx{H)^)ih d /A 3 , 
((E c L) 5 ( X (00) 52 ) 52 )i^/A 3 . (65) 
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For the leading order minimum given in Eq. ( 142|) . it is straightforward to see that the above 
higher dimensional operators only lead to corrections of 0(e 3 Vd) to m£ 22 . Hence, the NLO 
corrections to the charged lepton mass matrix dominantly arise from the modified vacuum 
alignment. Inserting the shifted values of x an d 4> given in Eq. (1621) into the leading order 
terms of we, we find that the charged lepton mass matrix can be parameterized as follows: 

2a 1 e 2 bis 2 b 2 e 

b\E 2 a 2 e —a\E 2 \ ev d , (66) 



mi 



a 2 e —a\E 
b 2 e —a\E 2 



in which the complex coefficients cjj and have 0(1) magnitudes. Consequently, the unitary 
matrix Ui, which is the transformation of the charged leptons to the physical basis in which 
m\mi is diagonal, is of the form 



U r 



1 

V- 

13 



and the charged lepton masses are 

1(2^ - b 2 2 /a 3 )e 3 v d \ 



iffy 



a 2 e v d \ 



\a>3£v d \ 



(67) 



(6* 



Hence, a realistic pattern of charged lepton masses can be produced. In the neutrino sector, 
the subleading operators are obtained by inserting the product X'P m au possible ways in 
the LO operators of w v and then extracting the A§ invariants. As we have shown above, 
the leading order structure of the neutrino flavons remains intact at NLO, and hence the 
corrections to the neutrino mass matrix only arise at NNLO. Therefore, the light neutrino 
mass matrix m u given in Eq. f )52|) and the diagonalization matrix U v of Eq. f )53|) are unchanged 
at NLO. Taking into account the corrections from the charged lepton sector, the lepton 
mixing angles are modified as follows: 

e 

i 



sin 9 



13 



(- 


_ M 


\a 2 





sin' 9 



12 



10 



sin 2 fl 



23 



1 

4 



h 

a 2 



03 



h 

a 2 

2 



«3 



a 1 
«3 



c.c. e 



+ c.c. 



(69) 



A^, while the deviation of the 



We see that #13 and 9\ 2 both receive corrections of 0(e) 
atmospheric angle #23 from its GR1 mixing value of ir/4 is only generated at 0(e 2 ). 

To investigate the predictions of our model in light of the current experimental ranges of 
the lepton mixing angles as described in the introduction and to cross-check the reliability of 
our analytical estimates, we have performed a numerical analysis in which the LO and NLO 
coefficients are taken to be random complex numbers with absolute value in the range of 
[1/3, 3], the corresponding phases are varied between and 2ti, and e is set to 0.04 (~ A^). 
The results of our analysis are shown in Fig{TJ We see that a large set of points fall into the 
region in which sin 2 #13 and sin 2 #12 are in the 3a interval, while #23 is in the la range for 
almost all the points. However, it is clear from the parametric analysis given above that a 
value of 6*i3 near the present upper bound (i.e., ~ A c ) would be unnatural in this model. 
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Figure 1: sin 2 Q\i and sin 2 #23 as a function of sin 2 #13. The horizontal and vertical lines represent the la, 
2cr and 3<r bounds on the mixing angles taken from Ref. [3]. 



5 Summary and Conclusions 

In this paper, we have studied the symmetry properties of lepton flavor mixing models in 
which at leading order the atmospheric mixing angle is 7r/4, the reactor angle #13 is zero, and 
the solar mixing angle is related to the golden ratio, (j) g — (1 + a/5)/2. We considered both 
of the proposed relations available in the literature: GR1 mixing, in which cot 6*12 = <fr g , and 
GR2 mixing, in which cos 612 = 4> g /2- We first investigated the symmetry properties of the 
neutrino mass matrices as given in a basis in which the charged leptons are diagonal, and 
found that in each case there is a Klein four (K4 ~ Z 2 x Z 2 ) symmetry that is preserved. 
Next, we investigated the question of determining the minimal (i.e., the smallest) discrete 
family symmetry group that encodes the symmetries of the fermion mass matrices. We found 
that A$ is the minimal group for GR1 mixing, and that up to a very high group order, there 
is no minimal discrete group that contains all of the symmetries needed for GR2 mixing. If 
the A 5 family symmetry is broken to Z 5 and in the charged lepton and neutrino sectors 
respectively, the GR1 mixing pattern follows immediately without fine-tuning. 

We then focused our attention on models of GR1 mixing based on A§ with this pattern 
of flavor symmetry breaking, and constructed two prototype supersymmetric lepton flavor 
models in which the SM leptons are all assigned to triplet representations of A 5 and GR1 
mixing is exactly reproduced at leading order. In the first model (Model 1), the complete 
family symmetry is A 5 x Z 5 , which is broken to Z 5 and as described. The model predicts 
the neutrino mass sum rule m]" 1 +m 2 _1 = 7713 , and the lightest neutrino mass is predicted to 
be larger than 0.011 eV (0.028 eV) for the normal (inverted) mass hierarchy. Furthermore, 
the Majorana phases and \m ee \ can also be expressed in terms of the light neutrino masses. 
However, although GR1 mixing is achieved without fine-tuning, the model has a drawback 
in that a separate fine-tuning of the charged lepton mass parameters is needed at leading 
order to achieve a phenomenologically acceptable mass hierarchy. 

To improve upon this situation, we constructed a second model (Model 2), in which 
the family symmetry is A§ x Z3 x Z%. In the neutrino sector, A§ is again broken to K4, 
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and both the normal and inverted hierarchy patterns of the light neutrino masses can be 
accommodated. In the charged lepton sector, the Z 5 is further broken, and the charged 
lepton mass hierarchies are controlled by the symmetry breaking parameters e ~ \ 2 such 
that there is no need to invoke a Froggatt-Nielsen U(l) symmetry. In this model, the tau 
and muon masses are of 0(evd) and 0(e 2 Vd) respectively, and electron is predicted to be 
massless at LO. We next studied the NLO corrections due to higher dimensional operators 
are in detail, and fount that the neutrino sector remains intact to NLO, but the charged 
lepton mass matrix receives corrections that generate small off-diagonal elements as well as 
the electron mass. To obtain an acceptable value for the electron mass, it is necessary to 
assume that the shift in the relevant flavon field value Sv^/v^ is suppressed at the level of 
0(e 2 ) instead of 0(e). The lepton mixing angles 9 12 and #13 are shifted from their GR1 
values by 0(e) ~ X 2 C while #23 receives corrections of 0(e 2 ) ~ A^. 

This model is consistent with current lepton data, as we verified explicitly by a numerical 
analysis. However, if the recent experimental hints of a potentially much larger value of #13 
are in fact verified in future experiments, this model would be strongly disfavored. Indeed, 
this tension with a larger value of 6*13 is a common feature within many lepton flavor models 
based on discrete family symmetries in which the reactor angle is zero at leading order (such 
as tri-bimaximal mixing models and others), as often the corrections should at most be of 
0(\ 2 ) to keep the solar mixing angle in a phenomenologically acceptable range. Hence, a 
precise measurement of #13 is a critically important test not only of this model, but also of 
many models within the general framework of discrete non-Abelian family symmetries. 
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Appendix: A§ Representations and Clebsch-Gordan Co- 
efficients 

As discussed in Section 2, we wish to work in a basis that satisfies Eq. ( fl5l) . i.e., one in 
which the order 5 generator T is diagonal. To determine this basis, we begin by noting that 
can be generated by two elements X and Y that satisfy the following relations [27] : 

X 2 = Y 5 = (Y 2 XY z XY' 1 XYXY~ 1 f = 1. (70) 

The explicit forms of these elements can be derived by carefully selecting several of the 
group representation matrices in the Shirai basis jl6], which has been done in [27] . For 
completeness, we list the explicit representations below, beginning with the the generators 
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of the three-dimensional irreducible representations. For the 3, X3 and Y3 are given by 



X, 



( 



\ 



4r \ 



4- 1 



Y« 



( 1 



V 



± 1 



(71) 



where as usual 4> g = (1 + \/5)/2 is the golden ratio. For the 3', the generators take the forms: 



A3' — — 



/ 



V 



1>g 
-1 



1 \ 



/ 



/ 



Y* = - 



1 \ 



For the 4-dimensional representation, we have 



X 4 



( -1 -1 -3 -V5 \ 

-1 3 1 -V5 

-3 1 -1 v% 

V -V5 -V5 V5 -1 / 



1 
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/ "I 
-1 

3 

V V5 
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V5 



-3 
1 
1 

V5 



v/5 \ 
-1 / 



and for the 5, the order two generator is 



X, 



4 
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2 

\ 40 s 



and the order five generator is 



Y K 
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fi 
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40 9 
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2 
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(72) 



(73) 



(74) 



(75) 



We note that Eqs. (174"]) - (175]) correct minor typos found in 

With the preceding representation matrices, it is straightforward to derive the following 
relationship between the Shirai basis and the generators S and T in Eq. ()15|) : 



S = XY 3 XY 2 X, T = Y 4 



(76) 



Applying the above transformation to the Shirai basis given above results in a T that is not 
diagonal but an S that is. Since we wish to work in the diagonal T basis, we diagonalize T 
via a unitary transformation U, and hence 



S = U ] XY' i XY 2 XU, T = U^Y 4 U, 



(77) 
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For the 3, this results in 
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with 



y/2 5 1 '* 
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• 5 l/4 _^l/4 
1 1 



(79) 
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where p 



e s 



The generators of the K± group given in Eq. (jSJ) can be expressed in terms 



of £3 and T 3 as follows: G\ = S 3 , Gi 
we have 
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and 
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V25 1 / 4 
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T^i S 3 iTq S 3 iT^i and G3 



l 3 ,D 3 'l 3 ,D 3 'l 3 ,D 3 i. 



This also contains K 4 , in which G 2 = S 3 r, G% 
Therefore, both the 3 and 3' irreducible representations in this basis for A 5 contain the 
matrices which generate the full set of symmetries of the neutrino mass matrix in the case 
of GR1 mixing. Turning now to the 4, we have 
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Finally, in the case of the 5, 
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and 
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1& /A ^g 
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(85) 



From the above representation matrices, it is straightforward to calculate the Clebsch-Gordan 
coefficients for the decomposition of the product representations, which we now report for 
this basis in detail. We will use oti to denote the elements of the first representation and $ 
to indicate those of the second representation of the product, and we use the subscripts "A" 
and "5"' to indicate a representation which is antisymmetric or symmetric, respectively. 
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V2a 1 f3 1 + a 2 f3 4 - a 3 (3 3 
V2oiif3 4 - a 2 (3 2 + a 3 (3i 

( + V2a 3 f3 3 \ 

a 1 f3 2 - V2a 3 (3 4 
-a 1 f3 3 + v / 2« 2 /?i 
\ -a x p A - V2a 2 (3 2 J 

V6a 2 (3 3 - VQa 3 p 2 \ 
\j2a x fi\ - 3a 2 (3 A - a 3 (3 3 
2V2 ai f3 2 + 2a 3 /3 4 

-2>/2ai/3 3 - 2a 2 /?! 
-V^ai/3 4 + a 2 /3 2 + 3a 3 /3i / 



3©5 = 3©3' ©4©5 
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-2«i/3i + y/3a 2 /3 5 + y/3a 3 /3 2 
y/3ai/3 2 + a 2 /3i - \/6a 3 /3 3 
y/3ai/3 5 - \ZQa 2 /3 4 + a s /3i 

/ V^cniPi + "2/^5 + a 3 (3 2 
3' ~ an/33 - y/2a 2 /3 2 - \/2a 3 (3 4 
\ ai/3 4 - V2a 2 (3 3 - V2a 3 (3 5 

( 2V2ai/3 2 - VQa^ + a 3 /3 3 \ 
—y/2a\fc + 2a 2 /3 2 - 3a 3 /3 4 
+ 3a 2 /3 3 - 2a 3 /3 5 

V -2 v / 2a 1 /3 5 - « 2 /3 4 + V^gfr / 

/ V^O^S - V^O^ \ 

-«i/3 2 - a/3« 2 /3i - a/2« 3 /3 3 
-2a 1 /3 3 - V2a 2 p 2 
2a^ 4 + v^a&ft 

V wife + y/2a 2 f3 A + v^a^i / 



3' 



y/Zaifii + a 2 /3 4 + a 3 /3 3 
ctxfc - y/2a 2 /3 5 - V2a 3 j3 4 
ai/3 5 - V2a 2 (3 3 - V2a 3 (3 2 

-2ol\P>\ + y/3a 2 /3 4 + y/3a 3 /3 3 
y/Sanfc + a 2 /3i - y/6a 3 /3 5 
\/3ai/34 - y/6a 2 /3 2 + a 3 (5i 

( V2~ ai f3 2 + 3a 2 p 5 - 2a 3 (5 4 \ 
2y/2a 1 f3 3 - y/6a 2 f3 1 + a 3 (3 5 
-2\p2a\fi 4 - a 2 /3 2 + ^/Qa 3 (5i 
- V2ai/? 5 + 2a 2 f3 3 - 3a 3 /3 2 J 

V3a 2 (3 4 - V3a 3 (3 3 \ 
2«i/3 2 + V2a 3 (3 A 
-aiP 3 - a/3« 2 /3i - \/2a 3 (3 5 
«i/3 4 + \/2a 2 (5 2 + y/3a 3 /3i 
-2ai/3 5 - V2a 2 (3 3 



3' 



I5 ~ «i/3 4 + a 2 /3 3 + a 3 /3 2 + a,$\ 

-a 1 (3 4 + a 2 (3 3 - a 3 (3 2 + a 4 /3i 
V2a 2 (3 4 - V2a 4 (3 2 
\f2onfo - V2a 3 ^ 
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3' A ~ 



«i/3 4 + a 2 (3 3 
s/2a 3 /3 4 
\/2«i/3 2 



«3/3 2 - a 4j 9i 



( a 2 f3 4 + a 3 (5 3 + a 4 /3 2 \ 
ai/?i + a 3 /3 4 + a 4 /3 3 
a 1^2 + a 2 /3i + a 4 /3 4 
\ ai/3 3 + a 2 /3 2 + a 3 /?i J 

I V3ai(3 4 - y/3a 2 (3 3 - V3a 3 /3 2 + a/3« 4 /3i ^ 
-V2a 2 (3 4 + 2V2a 3 (3 3 - V2a 4 (3 2 
-2V2ai/3 1 + \/2a 3/ 5 4 + V^a^s 
\/2ai/32 + v / 2« 2 /?i - 2^04/34 



3 ©3' 



3' 



2V2aiAs - \/2a 2/ 5 4 + \/2a 3 /3 3 - 2y/2atf 2 
—y/EaiPi + 2a 2 /3 5 + 3a 3 /3 4 - a 4/ S 3 
«i/3 4 - 3a 2 /3 3 - 2a 3 /3 2 + V^^Pi 

V^atfc + 2V2a 2 /3 4 - 2V2a 3 /3 3 - \/2a!4/3 2 
3ai/3 2 - y/§a 2 (3i - a 3 /3 5 + 2a 4 /3 4 
-2«i/3 3 + a 2 /3 2 + \/6o: 3 /3i - 3a 4 /3 5 

/ v^Mi - v^a 2 /3s + V^a 3 /?4 - 2v^a4/9 3 \ 
-V2ai/3 2 - \/3o; 2 /3i + 2V2a 3 (3 5 + V2a 4 (3 4 
V2 ai /3 3 + 2V2a 2 f3 2 - V3a 3 A - V2a 4 /3 5 
-2\/2ai/3 4 + v^a 2 /0 3 - V2a 3 (3 2 + x/3a 4 /?i / 



V 



5i 



/ y^fo - \/2a 2 A - V2a 3 (3 3 + y/2a 4 /3 2 \ 
—\[2ol\P\ — v3a 3 /3 4 — v^^^ 
A/3ai/3 2 + \/2a 2 f3i + \/3a3/35 
A/3a 2 /3 2 + A/2a 3 /3i + a/3« 4 /3 5 
— v3q!i/3 4 — A/3a 2 /33 — v2q! 4 /3i / 



V 



2ai/3 5 + 4a 2 /3 4 + 4a 3 /3 3 + 2a 4 /3 2 
4ai/3i + 2^2/35 
--y/6ai/3 2 + 2a 2 /3i - v^a^s + 2^/60!^ 
2a/6«i/33 - ^^2/32 + 2a 3 /3! - y/6a A (3 5 
2V6a 3 (3 2 + 4a 4 /3i 



\ 



5 g> 5 



>3', 



)4 A 



l s ~ ai/3i + a 2 /3 5 + a 3 /3 4 + a 4 /3 3 + a 5 /3 2 
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_ a 2 /3 5 + 2a 3 (3 4 - 2a 4 (3 3 - a 5 (3 2 
-\/3ai/?2 + V^a 2 f3 1 + A/2a 3 /3 5 - v^asAj 
y/3ai/3 5 + a/2«2/34 - A/2a 4 /3 2 - V^oi 5 /3i 



2a 2 /3 5 - a 3 (3 4 + a 4 /3 3 - 2a 5 /3 2 
v3ai/?3 — \/3a 3 [3i + V2~ot 4 /3 5 — y/2a 5 f) 4 
-v^ai/34 + a/2« 2 /3 3 - a/2« 3 /3 2 + y/3a 4 /3 



I 3V2ai/3 2 + 3a/2<*2/3i - V^asAs + 4^4^ - \/3a 5 /? 3 \ 
Sv^ai/^ + 4v / 3a 2 /3 2 + 3^(^/31 - y/3a 4 /3 5 - \/3a 5 f3 4 
3^01/34 - v / 3«2/?3 - v^a^ + 3a/2q;4/3i + 4^05/35 

V 3^1/^5 - \/3a 2 /? 4 + 4V3a 3 /3 3 - v^ft + Sy^ag/?! / 

/ x/2 ai /3 2 - V2a 2 /3 X + V3a 3 /3 5 - V3a 5 /3 3 \ 
^ ^ —V2on/3 3 + y/2a 3 /3 4 + y/3a 4 /3 5 - y/3a 5 /3 4 
-\[2axP> 4 - y/3a 2 /3 3 + y/3a 3 /3 2 + \f2a 4 (3i 
\ v^ai/Js - y/3a 2 (3 4 + A/3a 4 /3 2 - \f2a^ x / 

/ 2«i/3i + a 2 /3 5 - 2a 3 f3 4 - 2a 4 /3 3 + a 5 (3 2 \ 
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